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Abstract 

Tomography  using  CT  scans  and  MRI  scans  is  now  well-known  as 
a  medical  diagnostic  tool  which  allows  for  detection  of  tumors  and  other 
abnormalities  in  a  noninvasive  way,  providing  very  detailed  images  of  the 
inside  of  the  body  using  low  dosage  X-rays  and  magnetic  fields.  They 
have  both  also  been  used  for  determination  of  material  defects  in  mod¬ 
erate  size  objects.  In  medical  and  other  applications  they  complement 
conventional  tomography.  There  are  many  situations  where  one  wants 
to  monitor  the  electrical  conductivity  of  different  portions  of  an  object, 
for  instance,  to  find  out  whether  a  metal  object,  possibly  large,  has  in¬ 
visible  cracks.  This  kind  of  tomography,  usually  called  Electrical  Im¬ 
pedance  Tomography  or  EIT,  has  also  medical  applications  like  monitor¬ 
ing  of  blood  flow.  While  CT  and  MRI  are  related  to  Euclidean  geometry, 

EIT  is  closely  related  to  hyperbolic  geometry.  A  question  that  has  arisen 
in  the  recent  past  is  whether  there  is  similar  “tomographic”  method  to 
monitor  the  “health”  of  networks.  Our  objective  is  to  explain  how  EIT 
ideas  can  in  fact  effectively  be  used  in  this  context. 

1  Introduction  and  preliminaries 

Networks  have  become  ubiquitous  in  present  society  and  thus  it  has  become  important 
to  avoid  and  detect  disruptions.  In  particular,  it  is  important  to  prevent  malicious  in¬ 
truders  from  disrupting  them.  To  achieve  this  sufficiently  early,  it  is  essential  to  count 
on  a  mathematical  model  that  can  allow  early  detection  of  attacks  to  the  network.  The 
mathematical  tool  that  we  consider  to  accomplish  the  early  detection  of  disruptions  is 
based  on  the  use  of  tomographic  ideas.  One  of  the  questions  we  are  considering  is  how 
to  find  out  whether  an  attack  against  the  network  by  traffic  overload  is  taking  place  by 
monitoring  traffic  only  at  the  periphery  of  the  network  (input-output  map),  and  hence, 
the  use  of  a  tomographic  approach. 

1  This  write-up  is  an  extended  version  of  the  invited  lecture  given  by  the  third  author  at  PASI  2003.  It 
represents  research  supported  in  part  by  grants  N&F-DMS0070044  and  ARO-DAAD 19-01- 1-0494. 


We  present  first  a  general  idea  on  how  tomography  can  be  used  to  do  such  monitor¬ 
ing  and  implemented  as  a  diagnostic  tool.  Consider  subsets  A  and  Cl  in  R2  such  that 
A  C  Cl  and  let  /  €  C0(K2)  be  such  that  SuppfcCl.  Assume  A  represents  an  object 
from  which  we  want  to  get  some  information.  We  could  just  ask  for  /  to  be  rapidly 
decaying  for  2— >  00. 

Let  <1>  be  the  collection  of  all  the  straight  lines  in  Cl  connecting  any  pair  of  points 
a  and  b  where  a.b^dCl,  aj^b,  i.e.,  $  ={£  line  /a,b  €  £  fl  OCl,  a  ^  b },  and  /( x) 
represents  the  distribution  or  density  at  the  point  x  £  Cl.  Let  the  Radon  transform 


RU)  = 


f(x)dx 


(1.1) 


be  the  set  of  all  line  integrals  of  the  function  /. 

The  concept  of  tomography  then  can  be  understood  as  the  reconstruction  of  the 
function  /  from  the  set  of  values  given  by  R(f).  To  recover  /  is  thus  the  same  as 
finding  the  inverse  of  the  operator  R  .  Hence,  we  are  able  to  know  the  value  of  /  at  any 
point  x  in  Cl  without  having  access  to  the  interior  of  Cl. 

A  well-known  example  of  tomography  is  transmission  CT  in  diagnostic  radiology, 
[31],  where  CT  stands  for  computerized  tomography.  Essentially  the  setup  consists  of 
a  detector  and  an  X-rays  beam  source.  A  cross-section  of  the  human  body  is  scanned 
by  a  thin  X-rays  beam.  Because  the  density  of  the  tissue  of  the  human  body  changes 
from  its  surface  to  its  interior,  there  is  a  intensity  loss  which  is  recorded  by  the  detector 
and  processed  by  a  computer  to  produce  a  two-dimensional  image  which  in  turn  is 
displayed  on  a  screen.  Given  that  the  X-rays  go  through  the  tissue,  it  is  clear  the  X-rays 
absorption  is  related  to  the  attenuation  coefficient.  Let  f(x)  be  the  X-ray  attenuation 
coefficient  of  the  tissue  at  the  point  x.  Taking  a  close  look.  X-rays  traversing  a  small 
distance  Ax  at  x  suffer  a  relative  intensity  loss. 


A I 

~T 


fix)  Ax. 


(1.2) 


If  the  X-rays  are  considered  as  straight  lines,  as  indeed  they  essentially  are,  let  £  be 
the  straight  line  representing  the  beam,  ID  the  initial  intensity  of  the  beam,  and  1 1  its 
intensity  after  having  traversed  the  body.  It  follows  from  (1.1)  that 


j-  =  exp{-  J  f(x)dx}  (1.3) 

C 

thus  the  scanning  process  provides  us  with  the  line  integral  of  the  function  /  along 
each  of  the  lines  From  the  set  knowledge  of  all  of  these  integrals  the  problem  is  to 
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reconstruct  /.  Equally  well  known  by  now  is  MRI,  magnetic  resonance  tomography, 
where  the  underlying  space  is  ]R3  and  the  integrals  take  place  over  the  family  of  all 
planes  in  R3. 

See  [25]  for  other  examples  of  imaging  and  detecting  equipment  based  on  tomo¬ 
graphic  principles  and  [34]  as  a  recent  overview  of  this  kind  of  problems  discussed 
here.  We  refer  to  [5]  for  more  details  on  MRI. 

Sections  2  and  3  are  provide  the  background  information  on  the  Radon  transform. 
In  section  4,  we  also  consider  tomographic  examples  like  the  geodesic  Radon  transform 
in  the  hyperbolic  plane,  which  appears  naturally  in  relation  to  the  inverse  conductivity 
problem  and,  possibly,  to  internet  tomography.  We  conclude  in  section  5  with  some 
new  results  on  this  last  subject.  The  key  ingredient  is  the  attempt  to  understand  what 
happens  in  a  network  from  “boundary  measurements”,  that  is,  to  determine  whether 
all  the  nodes  and  routers  are  working  or  not  and  also  measure  congestion  in  the  links 
between  nodes  by  means  of  introducing  test  packets  (ICMP  packets)  in  the  “external” 
nodes,  that  is,  the  routers.  The  question  of  finding  out  whether  there  are  nodes  that 
are  in  working  order  is  a  classical  question  in  graph  theory.  For  networks,  it  is  also 
interesting  to  try  to  predict  future  problems  due  to  congestion.  (Note  that  nodes  could 
fail  to  work  for  other  reasons  than  congestion  on  the  links  starting  at  a  given  node.) 
This  requires  to  monitor  also  traffic  intensity,  also  known  as  load,  congestion,  etc.,  in 
different  contexts.  There  is  another  analogy  to  mathematical  tomography  that  arose 
independently  and  maybe  closer  to  the  consideration  of  this  question  in  the  context 
of  electrical  networks.  Curtis  and  Morrow  have  done  very  interesting  work  in  this 
context,  both  theoretical  and  in  simulations,  see,  for  instance,  [20]  and  [19].  Another 
analogy  in  the  same  direction  arises  when  we  consider  very  large  networks,  as  the 
internet,  which  could  be  considered  as  the  discretization  of  an  underlying  continuous 
model.  In  this  way,  we  can  see  the  analogy  with  the  well-known  inverse  conductivity 
problem  and  we  could  try  to  profit  from  the  large  body  of  mathematical  research  in 
this  area.  The  analogy  with  this  particular  inverse  problem  indicates  that  if  one  were 
to  pursue  this  “abstract”  approach  the  “correct”  geometry  is  closer  to  be  hyperbolic 
than  to  be  Euclidean  [7],  On  the  other  hand,  as  of  this  moment,  we  have  found  that 
those  tomographic  analogies  are  more  useful  for  providing  directions  of  research  and 
methods  to  consider  these  problems  than  providing  an  exact  correspondence  between 
the  two  phenomena.  It  is  in  this  context  that  [9]  modelled  “internet  tomography”  as 
an  inverse  Dirichlet-to-Neumann  problem  for  a  graph  with  weights.  In  this  situation, 
one  can  prove  that  characteristics  of  the  graph,  namely,  its  connectivity  and  the  traffic 
along  links  can  be  uniquely  determined  by  boundary-value  measurements  as  shown  in 
[9]  which  is  the  natural  analogue  of  the  continuous  inverse  conductivity  problem. 

Among  the  questions  that  arise  naturally  using  the  inverse  conductivity  problem  as 
a  guiding  model  there  are  a  number  of  questions  that  have  been  previously  addressed 
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using  other  points  of  view.  Namely,  the  problems  already  addressed  in  [17]  for  internet 
tomography  are: 

1 .  Link-level  inference,  in  other  words  link-level  parameter  estimation  based  on 
end-to-end  path-level  traffic  measurements.  Examples  of  this  are  unicast  inference  of 
link  loss  rates,  unicast  inference  of  link  delay  distributions,  topology  identification,  loss 
rates  by  using  multicast  probing  and  so  on. 

2.  Path-level  inference  (origin-destination  tomography  OD)  in  other  words  sender- 
receiver  path-level  traffic  intensity  estimation  based  on  link-level  traffic  measurements. 
One  example  of  this  is  time-varying  OD  traffic  matrix  estimation. 

We  would  like  to  conclude  by  thanking  the  editors  and  the  referee  for  his  useful 
comments. 

2  The  Radon  transform  in  M2 

Let  uj  £  S1,  then  u>  =  (cos  9,  sin  9),  and  take  pel.  The  locus  of  equation  x  ■  ui  =  p 
represents  the  line  l  that  is  perpendicular  to  the  line  r  passing  through  the  origin  and 
forming  an  angle  9  with  the  real  line  R.  If  B  is  the  intersection  of  l  and  r,  the  euclidean 
distance  d  (signed)  from  B  =  poj  to  the  origin  is  equal  to  p. 

Consider  a  nice  function  /  or  any  reasonable  function  /  ,  for  instance  /continuos 
and  compactly  supported,  then  consider  the  line  integral  with  respect  to  Euclidean  arc 
length  ds, 

OO 

Rf(u,p):=  J  f{x)ds=  J  f{x0  +  tu>±)dt.  (2.1) 

X-UJ—p  —  OO 

where  xa  is  a  fixed  point  in  l,  i.e.  it  satisfies  x0  •  ui  =  p,  and  uR  =  (cos  9,  —  sin  ff),  the 
rotation  of  c o  by  7t/2.  When  p  and  u>  range  over  R  and  S  Respectively,  we  get  all  of  the 
lines  in  R2.  Usually  xQ  is  taken  as  poj. 

The  map  /  — >  Rf  is  called  the  Radon  transform  and  Rf  is  called  the  Radon  Trans¬ 
form  of  /.  We  refer  to  [31],  and  others  there  in,  for  a  detailed  exposition  of  the  Radon 
transform.  Clearly  Rf  is  a  function  defined  oil  S 1  x  R,  i.e.  the  family  of  all  lines  in  R2 
with  the  compatibility  condition  mentioned  in  [8]: 

(Rf)(-w,-p)=Rf(w,p)  (2.2) 

Given  that  l  doesn’t  change  when  u>  and  p  are  changed  to  Au>  and  A p,  A  /  0,  A  €  R, 
then  the  Radon  transform  can  be  extended  from  S 1  x  R  to  R2  x  R.  The  pair  (Acu,Ap) 
is  identified  with  (w,p),  and  the  extension  of  the  Radon  transform  satisfies 

Rp(Xu>,  Xp)  =  Rp,((jj,p),  (2.3) 
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therefore,  the  Radon  transform  can  be  extended  as  an  homogeneous  function  of  degree 
zero,  a  very  important  property  of  this  transform. 

Rf  can  also  be  defined  as 


Rf(u,p):=  [  f(x)6(p  —  uj  ■  x)dx  (2.4) 

J  R2 

with  5  the  1 -dimensional  delta,  which  allows  us  to  obtain  easily  the  properties  below  in 
the  natural  coordinate  space  of  lines.  In  particular  for  =  (w  1,1x2)  with  |cu|  =  1,  then 


jyfd  dRf(u,p) 

a  )  —  Ui 


and 


dxi 

dRf(u,p) 

du>i 


dp 


(2.5) 


d 


=  ~g^R(xi  /)(W>P) 


hence,  it  follows  that  if  Pm  (, x )  is  a  homogeneous  polynomial  with  constant  coefficient 
of  degree  m  and  |w|  =  1, 


R(Pm(dx)  f(u,p))  =  Pm(w)  •  (2-6) 

and 

Pm{du)(Rf{w,p))  =  (-l)m-^R(Pm(x)  f(u,p)) 

where 

d  d  d  d 

and  obviously  x  =  (xi,  X2). 

One  can  similarly  define  the  Radon  transform  in  K"  and  verify  that  the  properties 
(2.5)  and  (2.6)  extend  to  this  case.  In  particular  for  the  Laplacian  A  in  R.  n, 


fl(A /) 

where,  for  each  direction  u>  £  Sn~ 1 
dimension  1 .  Note  that  in  general 


d2Rf(Lu,p) 

dp2 


(2.7) 


the  right  hand  side  is  the  Laplace  operator  in 


n/ a  r  \  (  \  _  /  2  I  1  2a ^ 

R\Af)(u,p)  —  (oj1  +  ■  ■  ■  ■  +ixn) - ^2 - 

As  a  consequence,  if  the  function  /  depends  also  on  time,  and  □  n  represents  the  wave 
operator  in  n  dimensions  we  conclude  that 


RDnf  =  D1Rf. 
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therefore,  the  Radon  transform  in  n  dimensions  is  localizable  if  and  only  if  the  wave 
equation  is  localizable.  Fixing  ut  €  S"  j  one  can  express  this  identity  by  saying  that 

o  2 

the  Radon  transform  interwines  the  wave  operator  □  „  =  A—  jfaj  in  n— dimensions  with 

o2  o2 

the  wave  operator  Di  =  in  1— space  dimension.  It  follows  that  the  Radon 

transform  can  not  be  localized  in  even  dimensions  [10].  In  spite  of  this  observation  one 
can  obtain  an  almost  localization  of  the  Radon  transform  in  R  2.  The  key  elements  is 
the  use  of  wavelets  as  it  will  be  described  in  the  next  section.  Meanwhile,  for  the  sake 
of  completeness  we  remind  the  reader  of  the  standard  inversion  formula  for  the  Radon 
transform  in  R2.  It  depends  on  the  following  identity,  usually  called  the  Fourier  slice 
theorem.  Namely,  writing  the  Fourier  transform  F2(f)  of  a  nice  function  /  in  R2  in 
polar  coordinates  (s,  to)  we  have 


/  R2 


f{x)  e“ 


^  dx  = 


Rf{u,p)e  lsp  dp,  x  £ 


or,  in  a  more  concise  form. 


(2.8) 


F2(f)  =  F(Rf) 

where,  clearly,  F2  stands  for  the  2-dimensional  Fourier  transform  and  F  stands  for  the 
1 -dimensional  Fourier  transform  in  the  variable  p  which  provides  one  standard  inver¬ 
sion  formula  for  the  Radon  transform 


/  =  F^F(Rf)  (2.9) 

There  is  another  inversion  formula  that  has  a  number  of  advantages  for  us,  and  we 
proceed  to  explain  it  now.  To  simplify  we  work  in  X  =  §(R2),  the  Schwartz  space  of 
functions  /  and  Y  =  S(S'1  x  R)  the  Schwartz  space  of  functions  g.  Let  fit  fa  £  X 
and  g\t  g2  £  Y,  and  (/i.  fa) x  ■>  (di,  52) Y  t^le  inner  products  in  X  and  Y  respectively, 
then  because  of  the  linearity  of  the  operator  II,  we  write  the  equation  that  defines  R  * , 
the  adjoint  operator  of  R 

(- Rf,g)Y  =  {faR*g)x  (2-10) 

The  explicit  expression  for  R*g  is  given  by 

I  g(u>,  u>  ■  x)du>  =  R* g  (2.11) 

Js1 

R*  is  called  the  backprojection  operator. 

The  function  R*g  is  such  that  for  x  fixed 


R*g( x)  =  /  g(u,  u>  ■  x)dio 
Js 1 


is1 

is  the  integral  of  g  over  all  lines  passing  through  x. 
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In  order  to  get  a  formula  for  /  from  the  Radon  transform  value,  the  next  important 
property  of  the  backprojection  operator  holds. 

(R*g)*f  =  R*(g®Rf)  (2.12) 

where  ©  stands  for  the  convolution  with  respect  to  the  second  argument,  and  the  Radon 
inversion  formula  is  then  given  by. 

f  ei2™-^F2(R*Rf )(<;)&  =  f(x)  (2.13) 

Jr2  1 

Introducing  the  square  root  A  of  the  Laplacian  operator  A,  we  have 

A  (R*Rf)(x)  =  f(x)  (2.14) 

which  is  usually  called  the  backprojection  inversion  formula . 

3  Localization  of  the  Radon  transform 

As  explained  above,  we  can  not  in  general  reconstruct  the  function  /  in  a  disk  D(a,  r ) 
of  M2  using  only  lines  l  passing  through  D(a,r).  One  can  localize  up  to  a  baseline 
value  of  the  function  /,  that  is,  one  can  recover  for  a  disk  /)(«,  r)  by  using  only  the 
data  Rf(l)  for  passing  through  D(a,r  +  e),  for  arbitrary  e  >  0,  up  to  an  additive 
constant  [11]  and  1 12]  .  The  key  element  is  the  use  of  wavelets. 

Let  us  recall  the  basic  properties  of  the  continuous  wavelet  transform  (CWT)  and 
the  discrete  wavelet  transform  (DWT). 

Let  Rbe  and  ft  the  translation  of  /  by  b,  i.e.  fb(x)  =  f(x  —  b),  then 

7fc(e)  =  e-i2ri7(c) 

where  /  is  the  Fourier  transform  of  /.  Now  let  Daf  be  the  dilation  of  /  by  the  scaling 
factor  a  £  M,  a  >  0  where  Daf  is  defined  as  Daf(x )  =  -^/(f )  where  the  term  -A 
is  chosen  such  that  ||/||2  =  ||-Da/||2  >  i-e-  /  and  Daf  have  the  same  energy,  then  one 
has 

(Daf)  (?)  =  D1/af(g)  =  Vaf(a ?)  (3.1) 

As  pointed  out  in  [27],  equation  (3.1)  tells  us  that  the  Fourier  transform  (Daf)  (?)  is 
dilated  by  1/a ,  then  we  lose  in  the  ?  -  domain  (frequency)  what  we  gained  in  x— domain 
(time).  In  other  words,  there  is  a  trade-off  between  time  and  frequency  localization  if  ? 
and  x  stand  for  frequency  and  time  respectively. 
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3.1  Wavelets  as  a  tool 

In  what  follows  for  /  £  Li(R)  (or  f  £  S(R) )  we  denote 


/(?)  =  /  /(*)  e-  i2^xdx,  (3.2) 

J  —  OO 

the  usual  Fourier  transform  of  /.  Let  us  recall  from  [27]  the  definition  of  the  continuous 
wavelet  transform  (CWT)  associated  to  a  “mother”  wavelet  'T.  Namely,  following  [8], 
given  a  “mother”  wavelet  'f  £  L2(R)  H  Li(K)  and  /  €  Z/2(M).  we  define  the  wavelet 
transform  of  /  as 

/•°°  —  t  —  b  dtt 

Wyf(a,b):=  /(f)' L( - )^=  =<  f,Da^b(t)  >l2  (3.3) 

J-oo  a  Va 

b,a£R,a>  0,  where,  for  a  function  g  and  b  £  M  we  let  gb(t)  =  g(t  —  b). 

One  requires  that  the  “mother”  wavelet  T  to  be  oscillatory,  i.e.  'If( x)dx  =  0. 
In  fact,  one  assumes  stronger  condition 


dg  <  oo, 


(3.4) 


called  the  admissibility  condition.  The  admissibility  condition  is  satisfied  when  T  has 
several  vanishing  moments, i.e.,  for  0  <  k  <  s 


xk  \I i{x)dx  =  0 


The  functions  Da^b  are  called  the  wavelets 

The  function  /  can  be  reconstructed  from  its  wavelet  transform  by  means  of  the 
“ resolution  identity”  formula 


f  =  C j,1  /  <  /,  Da*b(t)  >L2  Da*b(t)dt 

J  —  oo  J — OO 

where  Cy  <  oo  since  T  £  Li(R).We  refer  to  [27]  for  the  general  theory  of  wavelets. 
Proposition  1  explains  how  to  use  wavelets  to  obtain  (almost)  localization. 


Proposition  1  [10]  Let  n  be  an  even  integer,  and  h  £  L2(R)  a  function  with  compact 

support  such  that  for  some  integer  m  >  0  h  is  n+m-1  times  differentiable  and  satisfies 

1.  7  f/i  (7)  £  Li(R)  n  L2(R)  forO<  j  <  m,  0<  k  <  m  +  n  —  1 
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2.  t  ih(t)dt  =  0  for  0<  j  <  m+  1,  i.e.,  h  has  m  +  1  vanishing  moments 
Then 

1 1~nh(t)  =  o(\t\~n~m+1)as  |f|  i — >  oo 

and 

tn+m- i/i-^e  l2(K) 

The  fact  that  I  1~nh(t)  =  o(|t|  n  m+1)  as  |f|  i — >  oo  tells  us  that  I  1~nh  decays 
as  \t\~('n+m~1\  and  therefore,  it  does  a  good  localization  job. 

For  practical  purposes,  the  continuos  wavelet  transform,  CWT,  is  discretized  and 
the  discrete  wavelet  transform,  DWT,  is  obtained.  In  order  to  discretize  it,  consider 
m,  n  €  Z  and  the  values  a,  b  that  appear  in  W^f(a,  b)  are  restricted  to  only  discrete 
values  a  =  a™,  b  =  nb0a a0  >  1,  bQ  >  1  fixed.  (The  fact  that  a0  >  1,  bQ  >  1  it 
really  does  not  matter  because  to,  n  can  be  negative).  The  discrete  wavelet  transform 
DWT  of  /  is  defined  as 

/OO 

f(tW(a~mt-nb0 )  (3.5) 

-OO 

where,  as  before,  it  holds  that  f  'T(f)  =  0,  and  the  wavelets  are  given  by 

5 rm,n(x)  =  a~m,2^(a~mx  -  nbo) 

=  -  nb0a 

hence  'Tm  ri  is  localized  around  nb0a™  in  time,  (3.5)  can  be  also  expressed  as  (/,  T '  m,n) 
which  are  called  the  wavelet  coefficients. 

It  is  important  to  point  out  that  in  the  discrete  case,  in  general,  there  does  not  exist 
a  resolution  of  the  identity  formula  to  recover  /,  so  the  recovering  of  /  must  be  done 
by  using  some  other  means  ,  for  instance  numerical  ones.  The  choice  of  the  wavelet  T 
is  essentially  only  restricted  by  the  requirement  that  the  admissibility  condition  holds, 

oo  2 

i.e.,  =  f_OQ  — j^j — dg  is  finite.  Following  [21],  the  discretization  is  only  restricted 

to  positive  values  of  a  then  the  admissibility  condition  becomes 


,oo  £(?) 

C *  =  /  ——d^ 

Jo  kl 


— — — dg  <  oo 

kl 


and  since  a,  b  will  be  discrete  values  only,  then  the  dilation  parameter  is  chosen  as  a™, 
to.  €  Z  and  a0  1  is  fixed  (usually  a0  >  1).  The  value  b0  is  also  fixed  and  it  is  chosen 
such  that  'If(;r  —  n  ba)  covers  the  whole  line. 
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Now,  for  reasonable  'f  and  suitable  a0 ,  b0,  there  exist  \E'TOi„  so  that  the  discrete 
wavelet  coefficients  (  /,  *1 'm,n)  characterize  completely  /  which  is  given  by 

f  =  * m,n )  ^ m,n 

m,n 

then  any  function  in  L 2(H)  can  be  written  as  a  superposition  of  the  wavelets  *1 1  m,n- 

3.2  Wavelets  and  the  Radon  transform 

Now  we  want  to  state  some  results  that  relate  wavelets  and  the  Radon  transform,  which 
is  of  interest  for  tomography,  [11],  [12], 


Proposition  2  Let  p  £  L2(®0  teal  valued,  even,  and  satisfying 


(3.6) 


where  p  stands  for  the  1  -dimensional  Fourier  transform  of  p.  define  the  radial  function 
in  R2  by 


W?)  =  ^ 

where  as  before,  is  the  2-dimensional  Fourier  transform,  then  'h  is  a  wavelet  for  n=2 
and  the  wavelet  transform  offis  such  that 

Wvf(a,  b)  =  aF1!2  [  (WpRu  f)(a,buj)duj 

Js 1 

where  R u  f  is  such  that  R u  f(p)  =  Rf(w,p). 


Proposition  3  Let  'f  be  a  separable  2-dimensional  wavelet,  i.e., 
^(a;)  =  \E,1(xi)\E,2(a;2),  a:  =  (£1,0:2) 


where  for  i=l,2 


^(7) 


<  Ci(l  +  I7I)  1  for  all  7  £  R.  Defining  the  family  of  the 


one-dimensional  functions  (pu]}UjgS1  by 


1 


Pu  (7)  =  2  I7I  5,1(7^i)^2(7W2),  W  =  (wi,w2) 


10 


i-e.  Pu  =  F1  1(pul).  Then  for  every  f  £  Li(R2)  D  L2(R2), 

{W-q,f)(a,x)=a~1^2  f  (Wp  Rujf)  (a,x  ■  aj)dui 

Js 1 

The  proposition  shows  that  the  wavelet  transform  of  a  function  f(x)  given  any  mother 
wavelet  and  at  any  scale  can  be  obtained  by  backprojecting  the  wavelets  transform  of 
the  Radon  transform  of  /  using  wavelets  that  vary  with  each  angle,  the  argument  of 
cu,but  which  are  admissible  for  each  angle,  i.e.  C\\,  <  oo. 

4  The  hyperbolic  Radon  transform  and  EIT 

In  this  section  we  discuss  the  Radon  transform  on  the  hyperbolic  plane,  state  some  for¬ 
mulae  analogous  to  the  ones  that  were  given  in  section  2  to  invert  the  Radon  transform. 
The  backprojection  inversion  formula  is  one  of  them,  and  later  we  will  see  how  the 
hyperbolic  Radon  transform  is  related  to  electric  impedance  tomography  (EIT). 

In  [6]  and  [7]  it  is  shown  that  the  hyperbolic  Radon  transform  is  involved  in  the 
problem  of  reconstructing  the  conductivity  distribution  on  a  plate  by  using  electrical 
impedance  tomography  EIT. 

4.1  The  hyperbolic  Radon  transform 

Let  D  be  the  unit  disk  of  the  complex  plane,  i.e.  D  =  {z  £  C/  \z\  <  1}.  In  D,  a 
Riemannian  structure  is  defined  through  the  hyperbolic  metric  of  arc-length  ds  given 
by 


ds  = 


2  |  dz 


(4.1) 


with  dz  the  Euclidean  distance  in 
z,  w  £  D  is  given  by 


(1-10 

“2,  and  the  hyperbolic  distance  between  two  points 


7/  \  •  7  (  \Z~W\  \ 

diz.w)  =  arcsin  h  - -r - - 

U1  -  M  )1/2(i  -  M  )1/2/ 

The  set  of  lines  that  are  diameters  of  D,  and  the  set  of  intersections  between  the  Euclid¬ 
ean  circles  and  D  such  that  the  resultant  lines  (intersections)  are  perpendicular  to  the 
boundary  dD  of  1)  are  the  geodesics  or  h-lines  for  the  metric  (4.1). 

If  2  £  D  is  expressed  in  polar  coordinates  by  ( w ,  r)  where  w  =  z/  \  z\,r  =  d(z,  0), 
then  the  metric  (4.1)  becomes 


ds2  =  dr2  +  sinh2  r  dw2 


11 


where  dw2  is  the  usual  metric  on  dD.  This  indicates  that  the  area  in  hyperbolic  geom¬ 
etry  is  exponential  on  the  radius  r.  Let  us  recall  that  if  E  is  a  set  contained  in  the 
hyperbolic  disk  D ,  then  the  hyperbolic  area  of  E,  h  —  area(E),  is  given  by 


h  —  area(E)  =  /  /  dxdy 


- dxdy ,  z  =  (x,  y) 


IE  [(1-M2)]2 

and  the  hyperbolic  length  of  any  curve  7  in  D,h  —  length (7),  is  given  by 

2\dz\ 


h  —  length (7)  =  f 


J 7  i-l  *l2 

In  terms  of  the  Euclidean  Laplacian  A,  the  Laplace-Beltrami  operator  A  u  in  polar 
coordinates  on  D  can  be  expressed  as 


A  H  = 


U  U 

+  coth  r—  +  smh 
orz  or 


& 2 

r  dw2 


(4.2) 


and  in  the  Euclidean  coordinates,  z  =  (x,  y),  (4.2)  becomes 


A  H  = 


(1  —  x2  —  y2)2  .  d2  d2 
4  dx 2  dy2 


Following  [8],  we  denote  the  Radon  transform  of  a  function  /  within  the  hyperbolic 
plane  by  Rnf  which  is  a  function  on  the  family  of  geodesics  in  D.  It  is  defined  as 
follows, 

Rf{ 7)  =  RhUi)  =  [  f{z)ds(z),  7  geodesic  in  D,  (4.3) 

J  7 

where  /  is  a  function  such  that  (4.3)  makes  sense.  For  instance,  if  /  is  a  function 
with  compact  support,  i.e.,  /  £  C0(D)  or  f  £  §(Z?),  the  Schwartz  space,  which  is  the 
space  of  functions  rapidly  decaying  as  \z\  — »  1.  In  fact,  /  has  to  decay  a  bit  faster  than 
e~  r  because  the  element  of  length  ds(z)  grows  as  er.  (As  ds 2  =  dr2  +  sinh2  r  dw2, 
ds  grows  as  er).  If  I  is  the  space  of  all  the  geodesics  in  D,  then  the  dual  (adjoint) 
transform  R*  (backprojection  operator)  is  given  by 


#*$(2)  =  ®{l)dhz(Pt) 

where  T z  is  the  set  of  all  the  geodesics  passing  through  z  and  d/i  z  is  the  normalized 
measure  of  . 

Any  geodesic  passing  through  the  point  z  £  D  depends  only  on  one  of  the  end 
points.  (The  other  end  point  determines  the  same  geodesic  through  z.)  Therefore,  T  z  is 
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completely  determined  by  dD  =  Sl  hence  T z  can  be  parametrized  by  w  £  S 1  and  <7/x, 
is  naturally  associated  to  A.  dw. 

Having  done  this  parametrization,  the  purpose  now  is  to  invert  the  operator  R  // .  In 
order  to  invert  Rh  one  can  proceed  in  the  spirit  of  Radon’s  inversion  formula.  See  [24]. 
Following  [10],  one  try  to  find  a  filtered  backprojection  type  formula  like  2.13.  Recall 
that  if  k  £  Ljoc([0, 00))  we  can  associate  a  radial  kernel  such  that  for  f  £  C0(D ),  the 
convolution  operator  with  respect  to  this  radial  kernel  k  is  defined  as 


k  *  f(z )  =  k  *h  f{z)  =  [  f(w)k(d(z,w))dm(w)  (4.4) 

J  D 

where  dm(w)  is  the  measure  for  the  hyperbolic  area  which  in  polar  coordinates  is  given 
by 

dm  =  sinh  rdrdw 

Recall  the  formula  for  R*R  can  be  written  as 


R*Rf  =  R  *  /‘ 

The  analogous  result  for  the  hyperbolic  Radon  transform  Rh  is  given  by 

R*HRHf  =  k  *  f,  where  k(t)  =  — — — 

7T  sinh  t 

In  [7]  and  [8]  it  is  shown  that  by  letting  f  S(t)  =  coth  t  —  1,  we  obtain 

1 


4-7T 


Ah  S  *h  R*hRh  =  I, 


(4.5) 


(4.6) 


the  analogue  to  the  backprojection  inversion  formula  given  before. 

The  Fourier  transform  in  the  hyperbolic  disk  I)  for  a  radial  function  k  is  given  by 


k(X)  =  2ir  /  fc(f)PiA„i/2(cosh  t)  sinh  tdt,  for  Agl 

Jo 

where  Pu(r)  is  the  Legendre  function  if  index  //.  If  m  is  another  radial  function  then 


(k  *  m)(A)  =  k{\)m{\) 

as  we  know  it,  [8].  It  follows  that  as  k( A)  ^  0  VA  €  1  then  the  operator  Rh,  which 
takes  /  to  k  *h  f ,  is  injective. 
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4.2  Electrical  impedance  tomography  (EIT) 

EIT  has  a  number  of  applications  to  medicine  and  non-destructive  evaluation.  For  in¬ 
stance,  to  determine  the  existence  and  lengths  of  internal  cracks  in  the  wings  of  an 
airplane.  These  applications  are  related  to  the  inverse  problem  which  is  formulated 
now. 

Let  D  the  unit  disk  in  M2  and  (3  an  strictly  positive  function  defined  on  D  which 
is  unknown  and  represents  the  conductivity  distribution  inside  the  disk.  When  currents 
are  introduced  at  the  boundary  dD,  let  T  be  a  given  integrable  function  representing 
such  currents  and  such  that  the  average  of  the  values  of  T  on  dD  is  zero 


rkris  =  0 


idD 


and  consider  the  boundary  problem  with  Neumann  conditions 


(  div((3grad  u)  =  0,  in  D 
l  on  dD 

where  T  is  given  and  n  is  the  outer  unit  normal  vector  on  1)1).  This  problem  has  a 
unique  solution  n  where  the  uniqueness  of  a  is  up  to  an  additive  constant.  The  function 
u  is  the  potential  distribution  on  D  so  grad  u  is  the  electrical  field.  The  variation  of 
u  on  dD  has  to  correspond  to  the  known  values  of  T  on  01).  then,  if  s  represents  the 
tangent  vector  to  dD.  it  follows  that  the  tangential  derivative  of  u.  'jp- .  depends  linearly 
on  T.  So,  for  T  given  and  3.  the  unknown  conductivity,  there  exists  a  solution  u.  This 
defines  a  mapping 


P,  * 


du 

~dPs 


where  0  is  the  only  remaining  function  to  be  found. 

Let  Ap  be 

.  T  du 

P''  *  Ids 

Ap  is  a  linear  operator  from  the  Sobolev  space  H  a  ( dD )  into  Ha  1  (dD)  ,  and  / 3  deter¬ 
mines  A,;. Given  that  (3  is  to  be  found,  then  we  consider  the  nonlinear  mapping 


P  — >  A/3  (4.8) 

Now  the  problem  (Calderon’s  problem)  consists  in  determining  /3  once  A  p  is  given.  In 
other  words,  the  problem  is  to  find  the  inverse  of  the  mapping  (4.8),  and  this  problem  is 
called  the  inverse  conductivity  problem. 

Some  questions  arise  here.  Is  the  mapping  (4.8)  injective?.  If  so,  how  can  the  inverse 
of  A p  be  found?.  The  injectivity  of  it  in  two  dimensions  was  proven  by  Nachman  [29] 
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and  for  dimensions  higher  than  two  by  Sylvester  and  Uhlmann  [34]  and  [33],  See  [4] 
and  references  therein  for  some  uniqueness  results.  For  the  linearized  problem,  the 
injectivity  was  proven  by  Calderon  [13].  What  we  explain  next  is  how  to  try  to  find  an 
approximate  inverse. 

4.3  The  approximate  solution  to  the  EIT  problem 

As  it  was  stated  before,  0  is  called  the  conductivity  distribution  (of  some  object,  for 
instance,  lungs  tissue),  and  1/(3  is  called  impedance,  hence  the  name  of  EIT.  The  value 
of  (3  corresponding  to  different  constituents  like  human  lungs  tissue,  blood  and  so  on 
are  already  known,  then  one  only  looks  for  a  profile  of  the  areas  occupied  by  them.  EIT 
can  measure  the  rate  of  pumping  of  the  heart.  In  fact,  there  is  already  a  patented  device 
based  on  EIT  that  measures  that  rate. 

In  the  case  of  the  determination  of  cracks,  if  /3  >  0,  0  can  be  assumed  as  known, 
then  the  curves  and  the  existence  of  them  are  to  be  determined. 

We  want  to  emphasize  that  {3  is  generally  not  a  constant  but  in  the  case  we  consider 
now  (3  is  close  to  being  a  constant  positive  value  (3  ry  Assume  that  (30  is  initially  known, 
and  what  we  want  to  know  is  how  much  it  deviates  from  (3  0  where  (30 .  To  simplify  we 
set  (3a  =  1  so  the  deviation  of  0  is  governed  by 


0=1  +  6(3 

where  \S0\  «  1,  and  5(3  is  a  function  depending  on  the  position.  If  5(3  =  0  at  some 
point  w  in  the  object  being  studied  D.  then  there  is  no  any  “abnormal”  situation  at  w. 
It  is  also  assumed  that  there  is  no  any  deviation  on  <) D .  i.e.  5(3  =  0.  If  LI  is  the  solution 
of  (4.7)  for  0  =  1,  i.e., 

(  div(grad  U)  =  0,  in  D 
\  on  dD 

and  since  div(grad  U )  =  AU  (the  Laplacian  of  U),  it  follows  that 


(  AU  =  0,  in  D 
l  fr  =  *,  on  dD 


(4.9) 


now  let  u  be  the  corresponding  solution  of  (4.7)  for  the  perturbed  conductivity  0  = 
1  +  50  ,  then  there  is  a  perturbation  5U,  hence  u  =  U  +  5U.  The  perturbation  5U 
satisfies 

(  A(5 U)  =  —  ( grad  50,  grad  U) ,  in  D 
I  ff  =  ~m*>  on  dD 

and  since  T  represents  the  input  of  the  currents  and  they  can  be  arbitrarily  chosen  with 
the  only  constraint 


'Tris  =  0 


dD 
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then  the  input  SI/ ,  can  be  well  approximated  by  linear  combination  of  dipoles  where  a 
dipole  at  a  point  w  £  d D  is  given  by  Sw  the  Dirac  delta  at  w.  It  follows  that 

the  problem  (4.9)  for  the  dipole  (input)  —n  -§^8W  at  w  becomes 


AC/,,,  =  0,  in  D 


—S 

ds°w, 


dU,„ 

dn 


d  D 


(4.10) 


and  the  solution  Uw  of  (4. 10)  has  level  curves  which  are  arcs  of  circles  that  pass  through 
w  and  are  perpendicular  to  01).  Therefore,  the  level  curves  of  Uw  are  exactly  the  geo¬ 
desics  given  by  the  hyperbolic  metric.  At  this  point,  the  hyperbolic  Radon  transform  is 
involved  in  the  problem  and  can  be  used  to  solve  it. 

In  [7]  is  shown  that  the  linearized  problem  can  in  fact  be  described  explicitly  in  the 
context  of  hyperbolic  geometry  using  R  n  and  a  radial  convolution  operator  with  kernel 
k.  Let  k  be  given  by 

cosh_2(f)  —  3cosh~4(f) 


k(t)  = 


8i r 


then,  as  the  boundary  data  function  /j  =  d ^ '  -1  defined  on  the  space  of  the  geodesics  in 
D,  the  relation  between  8(3  and  //  can  be  shown  to  be 

Rn{k  *h  5(3)  =  n 

and  because  of  the  backprojection  operator,  one  obtains 


R*HRH{k  *H  5(3)  =  R*hh 

hence 

-7-  Ah{S  *h  {R*hh))  =  k  *h  8/3  (4.11) 

47 r 

Computing  the  hyperbolic  Fourier  transform  of  k,  k,  which  can  be  done  exactly,  it 
can  be  seen  that  k( A)  0,  VA  £  R,  and  consequently,  the  convolution  operator  with 
kernel  or  symbol  k,  k*H  is  invertible.  Formula  (4.1 1)  requires  to  invert  the  convolution 
operator  of  symbol  k  to  compute  5f3.  Barber  and  Brown  [2]  proposed  an  approximate 
inversion  and  Santosa  and  Vogelius  [32]  shows  that  the  inversion  formula  suggested  by 
[2]  is  a  generalized  radon  transform. 

To  numerically  implement  the  reconstruction  of  8(3  it  is  necessary  to  invert  the  geo¬ 
desic  Radon  transform  and  perform  a  deconvolution.  The  difficulty  of  numerically  im¬ 
plementing  (4.1 1)  lies  in  the  fact  that  it  is  complicated  to  numerically  implement  a  two- 
dimensional  non-Euclidean  convolution  on  the  hyperbolic  space.  In  [26],  Lissianoi  and 
Ponomarev  focus  on  the  problem  of  numerically  inverting  the  geodesic  Radon  trans¬ 
form  by  developing  an  algorithm,  and  the  problem  regarding  the  deconvolution  is  also 
considered  there.  For  this  purpose,  they  consider  the  inversion  formula  (4.6)  and  use 
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it  to  derive  an  inversion  formula  for  the  geodesic  Radon  transform  that  it  is  more  suit¬ 
able  for  computations.  The  interesting  open  problem  here  is  to  be  able  to  define  a  class 
of  “discrete  hyperbolic  wavelets”  that  provides  the  localization  described  in  section 
3  for  the  Euclidean  Radon  transform  and  has  computation  properties  similar  to  those 
of  the  Euclidean  wavelets.  For  examples  of  discrete  hyperbolic  wavelets,  we  refer  to 
[22]and[23], 

5  A  similar  problem  on  networks 

To  conclude  this  report,  we  would  like  to  mention  a  question  raised  by  the  second  author 
in  his  PASI  2003  lecture  and  the  new  results  obtained  since  then.  The  question  we  refer 
to  is  what  is  now  being  called  Internet  Tomography  [17]  and  [1].  The  problem  is  to 
be  able  to  find  out  whether  a  network,  usually  a  communications  network,  is  suffering 
some  sort  of  breakdown.  By  that  we  mean  that  traffic  along  the  network  either  can  not 
reach  every  node  in  the  network,  or  when  we  add  a  measure  of  traffic  around  nodes,  the 
traffic  is  so  large  in  some  parts  of  the  network  that  it  would  take  very  long  to  go  from 
one  node  to  another.  When  the  network  is  large,  the  information  is  naturally  gathered 
at  the  “periphery”  of  the  network  and  hence  the  name  of  internet  tomography.  The 
similarity  to  usual  tomography  becomes  closer  when  one  uses  as  a  way  to  measure  the 
traffic  “packets”  sent  from  the  “boundary  of  the  network”  and  measures  whether  they 
arrive  to  the  other  boundary  points  and,  more  often,  how  long  it  takes  to  get  there. 

Computer  scientists  have  done  “experimental”  work  on  this  subject  and  have  sug¬ 
gested  that  the  natural  model  of  internet  tomography  is  a  graph  situated  in  a  portion  of 
the  2-sphere  or,  what  is  essentially  the  same  thing,  in  the  hyperbolic  plane  [28].  Before 
we  proceed  further,  let  us  note  that  we  have  alluded  to  two  natural  types  of  “disrup¬ 
tions”  of  the  network.  First,  when  thought  as  a  planar  graph,  if  a  node  or  collection  of 
nodes  have  ceased  to  exist  because  of  an  “intrusion”,  the  “topology”  of  the  network  has 
changed.  There  has  been  very  significant  work  on  this  direction  by  experts  on  graph 
theory.  The  important  work  of  Fan  Chung  and  her  collaborators  offers  crucial  insights 
into  this  question.  (See,  for  instance  [14],  [15]  and  [16].)  Another  situation,  that  re¬ 
sembles  more  what  “conventional”  tomography  is  supposed  to  help  with,  arises  when 
traffic  among  certain  nodes  starts  to  increase  to  levels  where  the  graph  structure  remains 
intact  but  there  is  significant  slow  down  due  to  this  large  amount  of  traffic.  Communica¬ 
tion  networks  and,  regretfully,  road  networks  are  a  well  known  example  of  this  second 
phenomenon.  In  either  case,  the  desire  is  to  be  able  to  detect  this  problem  when  it  is 
incipient  to  try  to  devise  a  solution  to  it.  It  is  the  latter  problem  that  is  of  interest  to  us. 
One  can  see  that  Munzner’s  suggestion  leads  to  a  question  closely  resembling  EIT,  and 
it  is  natural  to  consider  it  a  problem  in  hyperbolic  tomography  of  the  kind  described 
earlier.  On  the  other  hand,  we  have  just  obtained  a  significant  result  on  the  inversion  of 
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the  Neumann-Dirichlet  problem  by  studying  it  directly  on  “weighted”  graphs  [9]. 

Let  us  explain  now  a  bit  more  in  detail  what  these  recent  results  are  and  what  new 
questions  they  open  up.  To  understand  the  ideas  better  let  us  consider  a  very  simple 
example  of  a  planar  network,  the  square  network  G  [20] .  This  network  is  constructed 
as  follows.  The  nodes  of  G  are  the  integer  lattice  points  p  =  (i,j)  with  0  <  i  <  n  +  1 
and  0  <  j  <  n+1  and  exclude  the  points  (0,0),  (n+1, 0),  (0,  n+1),  and  (n+1,  n+1). 
Let  V  be  the  set  of  nodes,  and  int  V  the  interior  of  V  consisting  of  the  nodes  p  =  (i,  j) 
with  0  <  i  <  n  and  0  <  j  <  n.  The  boundary  of  G  is  denoted  di l  and  it  is  equal 
to  V  \  intV.  Let  p  he  a  node  then  it  has  four  neighboring  nodes  which  are  the  nodes 
at  unit  distance  from  p.  Call  the  set  of  these  neighboring  nodes  as  N(p).  If  p  is  an 
interior  node  then  N (p)  is  in  V,  and  if  p  is  on  <9f i  then  it  has  only  one  neighboring  node 
which  is  the  interior  node  that  has  unit  distance  from  p.  If  a  line  segment  l  connects 
a  pair  of  neighboring  nodes  p  and  q  in  intV  or  if  it  connects  a  boundary  node  p  to  its 
neighboring  interior  node  q  is  called  edge  or  conductor  and  denoted  pq.  In  the  case  in 
which  p  is  on  the  boundary,  the  edge  is  called  a  boundary  edge.  The  set  of  edges  is 
denoted  by  E,  and  usually  the  graph  G  is  denoted  by  G(V. ,  E). 

Let  to  a  non-negative  real-valued  function  on  E,  the  value  n(pq)  is  called  the  con¬ 
ductance  of  pq  and  1  /u>(pq),  the  resistance  of  pq,  and  in  is  the  conductivity  (in  is  also 
called  a  weight).  A  function  u  :  V  — >  M  gives  a  current  across  each  conductor  pq 
by  Ohm’s  law,  I  =  n(pq)(u(p)  —  u(q))  (I  the  current).  The  function  u  is  called  in¬ 
harmonic  if  for  each  interior  node  p, 

^2  w(m)(m(<?)  -  u(p))  =  0 

qeN(p) 

then  the  sum  of  the  currents  flowing  out  of  each  interior  node  is  zero,  and  this  is 
Kirkhoff’s  law.  Let  <!>  a  function  defined  at  the  boundary  nodes,  the  network  will  ac¬ 
quire  a  unique  n-harmonic  function  u  with  u(p)  =  ^(p)  for  each  p  £  dG  in  other 
words,  <1>  induces  u  and  u  is  called  the  potential  induced  by  <L  Considering  a  con¬ 
ductor  pq  then  the  potential  drop  across  this  conductor  is  A u(pq)  =  u(p)  —  u(q). 
The  potential  function  u  determines  a  current  l  <\>(p)  through  each  boundary  node  p,  by 
/<j> (p)  =  uj(pq)  (u(p)  —  u(q)  ,  q  being  the  interior  neighbor  of  p.  As  in  the  continu¬ 
ous  case,  for  each  conductivity  in  on  E,  the  linear  map  A  from  boundary  functions  to 
boundary  functions  is  defined  by  =  /,],  where  the  boundary  function  <I>  is  called 
Dirichlet  data,  the  boundary  current  is  called  Neumann  data,  and  the  map  Au  is 
called  the  Dirichlet-to-Neumann  map. 

The  problem  to  consider  is  to  recover  the  conductivity  n  from  A  :jJ,  which  is  analo¬ 
gous  to  the  the  inverse  problem  in  the  continuous  case.  The  two  basic  problems  are  the 
connectivity  and  conductivity  of  the  network.  Note  that  the  connectivity  of  the  network 
or  the  situation  where  the  network  remains  connected  but  some  edges  disappear  is  a 
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topological  problem,  the  configuration  of  the  graph  has  changed.  For  detailed  theory 
about  electrical  networks,  planar  graphs,  recovering  of  a  graph  and  harmonic  functions, 
we  refer  to  [18]  and  the  work  of  Curtis  and  Morrow  [19]. 

The  discrete  or  finite  nature  of  graphs  makes  working  on  graphs  basically  easier  than 
investigating  these  problems  in  the  continuous  case,  although  it  gives  rise  to  several 
disadvantages.  For  example,  solutions  of  the  Laplace  equation  for  graphs  have  neither 
the  local  uniqueness  property  nor  is  their  uniqueness  guaranteed  by  the  Cauchy  data, 
contrary  to  the  continuous  case  where  they  are  the  most  important  mathematical  tools 
used  to  study  the  inverse  conductivity  problem  and  related  problems  [9].  The  inverse 
problem  that  we  study  is  to  identify  the  connectivity  of  the  nodes  and  the  conductivity 
on  the  edges  between  each  adjacent  pair  of  nodes. 

Given  a  network  with  a  pattern  of  traffic  measured  as  the  “usual”  load  between 
adjacent  nodes  (e.g.,  number  of  messages)  one  can  associate  to  it  a  Laplace  operator 
denoted  Aw,  where  the  weight  u>  is  a  sequence  of  values  representing  the  usual  loads 
between  every  pair  of  adjacent  nodes  in  the  network. 

We  define  the  degree  dux  of  a  vertex  in  the  weighted  graph  G  with  weight  u>  by 

dux  =  ^2  u;(x,y) 

v&v 

and  the  weighted  u-Laplacian  A uf  by 


A  uf  ■=  ^2[f(y)  -  xgV 

yev  a“X 

A  graph  S  =  S(V' ,  E')  is  said  to  be  a  subgraph  of  G(E,  V)ifVt<zV  and  Et  C  E. 
In  this  case,  we  call  G  a  host  graph  of  S.  The  integration  of  a  function  /  :  G  — >  R  on  a 
graph  G  =  G{V,  E)  is  defined  by 

/  /  :=  y2  f(x)dux  or  simply  /  fd„ 

Jg  xev  Jg 

For  a  subgraph  S'  of  a  graph  G  =  G(V,  E)  the  (vertex)  boundary  dS  of  S  is  defined  to 
be  the  set  of  all  vertices  z  €  V  not  in  S  but  adjacent  to  some  vertex  in  S,  i.e.. 


dS  :=  {z  £  V  |  z  ~  y  for  some  y  G  S} 

O 

and  we  define  the  inner  boundary  dS  by 


dS  :=  {z  €  S  |  y  ~  z  for  some  y  G  dS} 
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Also,  by  S  we  denote  a  graph  whose  vertices  and  edges  are  in  S  U  <9, S'. The  (outward) 
normal  derivative  -§^{z)  at  z  £  dS  is  defined  to  be 

:=  'Em  - 

n  yes  z 

where  d’u  z  =  £ye5w(z,y) 

An  attack  by  saturation  corresponds  to  a  new  weight  ui '  where  the  load  on  each  edge 
has  either  remained  the  same  or  increased  (substantially  in  some  parts  of  the  network). 

Associated  to  the  weight  ui  there  is  a  Laplace  operator  in  the  network  we  consider 
the  response  to  diagnostic  “proves”  applied  to  the  outside  boundary.  The  boundary 
observations  (outputs)  correspond  to  the  Neumann-to-Dirichlet  map  for  the  Laplacian 
Aw. 

The  theorem  below  shows  that  the  Neumann-to-Dirichlet  map  for  A  is  different 
to  that  for  Aw. 


Theorem  4  [9]  Let  lo i  and  LO2  be  weights  with  <  UJ2  on  S  x  S,  G  a  graph  andfi, 

f2  :  S  — *  K  be  functions  satisfying  that  for  j=l,  2, 

(  A =0,  x  G  S 

l  3^00  =  *(*), 

1  Is  fjduj  =  A 

for  a  given  function  $  :  dS  — >  K.  with  J9S  $  =  0,  and  for  a  suitably  chosen  number 
K  >0. 

If  we  assume  that 

O 

(i)  u>i  (z,  y)  =  u>2 (z,  y)  on  dS  x  dS 
(H)fl\dS  =f2\dS,  then  we  have 


fi  =  on  S 


and 


u>i  =  u>2  on  S  x  S 

whenever  fi(x)  ±  fi(y)  and  f2(x)  f2(y)  ■ 


Note  that 

j  A uf(x)  =  0,  x  G  S 

I  =  $(*)>  z^ds 

is  known  as  the  Neumann  boundary  value  problem  NBVP.  In  [9]  it  is  shown  that  the 
NBVP  has  a  unique  solution  up  to  an  additive  constant. 
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The  second  conclusion  of  the  theorem  shows  not  only  whether  or  not  each  pair  of 
nodes  is  connected  by  a  link,  but  also  how  nice  the  link  is.  Moreover,  the  proof  gives 
an  algorithm  to  detect  if  the  weights  change  on  the  edges. 

The  conditions  u>i  <  u)2  (monotonicity  condition)  and  JG  fjdUj  =  K  (the  normal¬ 
ization  condition)  are  essential  for  the  uniqueness  of  the  result  as  the  following  example 
from  [9]  shows. 

We  know  that  the  NBVP  has  a  unique  solution  up  to  an  additive  constant.  There¬ 
fore,  the  Dirichlet  data  f  \gs,  z  £  OS  is  well-defined  up  to  an  additive  constant.  Here 
we  discuss  the  inverse  conductivity  problem  on  the  network  (graph)  S  with  nonempty 
boundary,  which  consists  in  recovering  the  conductivity  (connectivity  or  weight)  oj  of 
the  graph  by  using  the  Dirichlet-to-Neumann  map  with  one  boundary  measurement.  In 
order  to  deal  with  this  inverse  problem,  we  need  at  least  to  know  or  be  given  the  bound¬ 
ary  data  such  as  f(x),  -§^{z)  for  z  £  dS  and  u>  near  the  boundary.  So  it  is  natural 
to  assume  that  /  |  ,  -§^\  and  u>\  0  are  known  (given  or  measured).  But  even 

dS  x  dS 

though  we  are  given  all  these  data  on  the  boundary,  we  are  not  guaranteed,  in  general, 
to  be  able  to  identify  the  conductivity  u>  uniquely.  To  illustrate  this  we  consider  a  graph 
S  whose  vertices  are  { 1,2,3}  and  dS  =  {0, 4}  with  the  weight 

w(0, 1)  =  1,  w(0,  k)  =  0  (k  =  2, 3, 4), 

and 

w(3, 4)  =  1,  u(k, 4)  =  0  (k  =  0, 1, 2). 


Let  /  :  S  — >  R  be  satisfying  A uf(k)  =  0,  k  =  1,2,3.  Assume  that  /( 0)  =  0, 

/( 1)  =  1,  /( 3)  =  3,  /( 4)  =  4,  /( 2)  =  unknown.  Thus,  since  dS  =  {1,3},  the 

boundary  data  f  ,  4^-  a  and  cj  I  „  are  known  . 

J  •'WtLn'as  'asxa°s 


In  fact. 


the  problem  is  to  determine 


w(l,  2)  =  x,  u>( 2,  3)  =  y,  uj(  1, 3)  =  2,  and  /( 2). 

From  Aw/(fe)  =  0,  k  =  1, 2,  3,  we  have 

,m  _  f(0)+xf  (2)+3:  __  1 
■I  G?  l+i'  +  :  ’ 

=  xf(1')+yf(3') 

ft  o')  __  i/(l)+y/(2)  +  /(4)  _  q 
J  \  /  z+y+1 

This  system  is  equivalent  to 

tit  /  x(y-l)+y{x~l)  +  2z(x  +  y)  =  0, 

(1)  l  m  =  fg- 

This  systems  has  infinitely  many  solutions.  For  instance,  assume  that  z  =  0,  that  is,  the 
two  vertices  1  and  3  are  not  adjacent.  Then  (1)  is  reduced  to 

r  -  +  -  =  2 

(2)  \/(2)  =  i@t 

It  is  easy  to  see  that  there  are  infinitely  many  pairs  (x,  y)  of  nonnegative  numbers  sat¬ 
isfying  the  first  equation  in  (2),  so  that  /( 2)  is  undetermined  as  a  result.  In  view  of 
the  above  example,  in  order  to  determine  the  weight  u ;  uniquely  we  need  some  more 
information  than  just  f\dSi  -§^\dS  and  a>|  o  .  If  we  impose  in  this  example  the 
additional  constraints  that 

x  >  1,  y  >  1  and  n  >  0 

in  (1),  then  the  equation  fl)  yields  a  unique  triple  of  solution  x  =  1,  y  =  1,  z  =  0  and 

/( 2)  =  2. 

There  are  many  problems  to  be  answered,  for  instance  what  happens  if  the  number 
of  nodes  is  not  finite?  What  is  the  hyperbolic  version  of  the  discrete  case?.  If  we  allow 
to  consider  also  to  =  0  then  the  presence  of  zero  weights  tells  us  that  the  conductivity 
on  the  edge  ( a  particular  one)  is  either  down  or  the  nodes  connected  to  that  edge  “disap¬ 
pear”  in  the  sense  that  the  edge  length  becomes  infinite  and  this  is  because  uniqueness 
is  not  true.  We  still  need  to  get  stronger  results  to  determine  the  configuration  of  a 
network  (connectivity).  Let  us  add  that  very  recently  Bensoussan  and  Menaldi  [3]  have 
given  a  slightly  different  proof  of  theorem  4  relying  on  the  fact  that  A  u  is  a  positive 
operator. 
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